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ABSTRACT 

Context. We report on results of nearly two years of INTEGRAL/SPI monitoring of the Galactic microquasar GRS 1915+105. 
Aims. From September 2004 to May 2006, the source has been observed twenty times with long (~100 ks) exposures. We present an 
analysis of the SPl data and focus on the description of the high-energy (> 20 keV) output of the source. 
Methods. Temporal and spectral analysis of the SPl data. Comparison to simultaneous 1.2 ~ 12 keV ASM data. 
Results. We found that the 20-500 keV spectral emission of GRS 1915+105 was bound between two states. It seems that these 
high-energy states are not correlated with the temporal behavior of the source, suggesting that there is no direct link between the 
macroscopic characteristics of the coronal plasma and the the variability of the accretion flow. All spectra are well fitted by a thermal 
comptonization component plus an extra high-energy powerlaw. This confirms the presence of thermal and non-thermal electrons 
around the black hole. 
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1. Introduction 

GRS 1915+105 is among the most notorious accreting black 

■ holes in our Galaxy. Not only is it one of the brightest and 
most variable X-ray sources in the sky (see Castro-Tirado et 
al. 1992 for a first detection report and Belloni et al. 2000 for 
detailed variability analysis), but it is also the first Galactic 

■ object in which superluminal plasma ejections have been ob- 
served (Mirabel & Rodriguez 1994). Consequently, much work 
has been done to try to understand the complex nature of the 
accretion flow (see e.g. the comprehensive analysis by Done, 

. Wardziriski & Gierliriski 2004) and the disk-jet coupling in this 
' binary system. An extensive review on these aspects can be 
' found in Fender & Belloni (2004). Very recently, Rodriguez 
et al. (2008a & 2008b) presented an extended report on GRS 
1915+105 discussing insights gained from a 2-year long moni- 
toring campaign involving multi-wavelength observations. 

However, most of the recent studies focus on the astonish- 
ing X-ray properties, the soft y-ray emission (> 100 keV) being 
generally observed at poor signal to noise ratio. Understanding 
the high-energy behavior of the source is nevertheless very im- 
portant as it is assumed to trace the physical processes occurring 
at the innermost regions surrounding the black-hole (e.g. Galeev 
et al. 1979 or Malzac 2007 for a review). Featuring good spec- 
tral resolution and sensitivity up to several MeV, SPl is a good 
instrument to tackle this issue. 

The emission ranging from hundreds of eV up to almost 
1 MeV is generally assumed to originate from two components 
(e.g. ViUiu et al 2001). While the softer part is unanimously at- 
tributed to thermal emission from the disk, various interpreta- 
tions exist concerning the origin of the high-energy emission. 
Detailed analysis of joint RXTE/OSSE observations allowed 
Zdziarski et al. (2001 & 2005) to rule out some of these inter- 
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pretations. They suggested that the observed 10-500 A:eV emis- 
sion originated from inverse Compton scattering of the soft disk 
photons by a hybrid thermal/non-thermal electron population. 

In this paper we present all SPl observations of GRS 
1915+105 from September 2004 to May 2006. We focus on 
gaining the most accurate high-energy picture of the source, 
mainly through extensive spectral analysis. First we briefly de- 
scribe the instrument as well as our data reduction methods. The 
observational results are presented and further discussed in the 
concluding part of the paper. 

2. Observations and data reduction 

2.1. INTEGRAL/SPI 

SPl is a high resolution y-ray spectrometer (Vedrenne et al. 
2003) aboard the INTEGRAL observatory. The observational 
strategy of the INTEGRAL mission is based on approximately 
3-day long revolutions during which one or several fields of view 
are sampled by means of 30-40 minute long fixed pointings, 
separated by a 2° angular distance. As the number of SPl pixels 
is small (among the initial 19 detectors, two broke down reduc- 
ing the detector plane to only 17 pixels), this so-called 'dither- 
ing' scheme (see Jensen et al. 2003 for details) is essential for 
SPl image reconstruction. In particular, combining the data from 
a set of dithered pointings allows one to enhance the precision 
of flux extraction (see e.g. Joinet et al. 2005). 

2.2. Data analysis 

We analysed all public data from nearly two years of SPl ob- 
servations on GRS 1915+105 (September 2004 - May 2006). 
For each revolution, we first gathered all pointings where GRS 
1915+105 was less than 12° olf the central axis. Pointings show- 
ing contamination by solar flares or radiation belt exit/entry were 
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excluded. The log of the resulting 1 .7 Ms of observational cov- 
erage is given in Table 1 . In order to detect the emitting sources 
in the field of view, we used the SPIROS software (Skinner & 
Connell 2003) to produce 20-50A:ey images for each observa- 
tion. The positions of the active sources (detected at a minimum 
of 5 o") were then given as input to a specific flux-extraction al- 
gorithm, using the SPI instrument response for sky-model fit- 
ting. From there we started to analyse the background and source 
behavior for each observation. We first allowed the background 
pattern as well as the main source(s) to vary between successive 
pointings. In this way we determine each component's most ap- 
propriate variability timescale for the final data reduction. We 
used pointing-durations (x; 2ks) as a timescale for the GRS 
1915H-105 light-curves and adapted the time sampling for sub- 
sequent extraction of the 20 keV -8 MeV spectra according to 
the observed temporal behavior. This method allowed us to min- 
imize the error bars without losing any scientific information. 

3. Results 

3.1. Light-curves 

Figure 1 displays the total 20-50 keV light-curve for our obser- 
vational period. GRS 1915+105 shows relevant long-term vari- 
ability, the averaged source flux per observation one-day) 
spans between 90 and 3E0 mCrab, with an approximate uncer- 
tainty of 5 mCrab. On the shorter science-window timescale 
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Fig. 1. Long-term variability of GRS 1915+105 in the SPI band. 
Each point represents the mean 20-50 keV flux per observation 
(corresponding to a timescale of about one-day, see Table 1). 



{^Iks), the obtained individual light-curves show in general 
lower variability. Within a single observation, the source flux 
varies by at most a factor of 2. As GRS 1915+105 is well 
known for being very variable in X-rays, we also considered 
SPI-simultaneous 1 .2 - 12 keV ASM light-curves to compare the 
X- and soft y-ray behavior of the source. We found an anti- 
correlation trend between the ~ one-day averaged 1.2-12 and 
20 -5QkeV source fluxes (Figure 2) with a linear correlation fac- 
tor p of -0.59 + 0.02 and a 99% significance. On the science- 
window timescale however, a positive flux-correlation is dis- 
cernible for observation 368, as the source exhibits very strong 
variability (Figure 5 right). The variability amplitudes (param- 
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Fig. 2. One day timescale flux-flux relation between the 1 .2 - 12 
and 20 - 50 keV bands. We observe an anticorrelation in the flux 
trend between the two bands. 



eterized by fractional rms) are not straightforward to compare 
since the time bins are very different between the two instru- 
ments. More importantly, the Icr SPI errors for the individual 
light-curves are of the same order of magnitude than absolute 
rms variability, thus limiting a solid scientific interpretation. 

3.1.1. Spectral analysis 

We used XSPEC 11.3.2 (Arnaud 1996) for spectral analysis. 
First, each observation-averaged spectrum was fitted with a ba- 
sic powerlaw model which gives a good description of most of 
the data. Then we characterized the spectra by 20-50A:ey flux 
and best-fit photon index, the latter being found to range from 
2.8 to 3.5, with typical uncertainty of 0.1. The results are sum- 
marized in Table 1 and illustrated by a photon index versus flux 
diagram given in Figure 3. We used this figure to select four 
observations which will be presented in more detail in the fol- 
lowing subsections. We selected the observations made during 
INTEGRAL revolutions 295 and 423 since they both have rather 
high average fluxes (and therefore good signal to noise ratios), 
but at the same time very different spectral shapes. Further, the 
light-curves showed no significant variability during these obser- 
vations (neither in the 20 - 50 ^eV band nor in X-rays), hence al- 
lowing a meaningful analysis of the observation-averaged spec- 
trum. In addition, we selected two other observations (Rev. 246 
and 368) because they were part of large multi-wavelength cam- 
paigns, allowing us to put the SPI data in a broader context. 

3.1 .2. Observations 295 and 423 

From the ASM light-curves in Figure 4, we see that during both 
observations GRS 1915+105 showed very similar X-ray activity, 
characterized by low flux and almost no variability. The tempo- 
ral properties in the 20 -5QkeV band are also quite similar, with 
again little variability during both observations. On the other 
hand, spectral characteristics are found to be significantly differ- 
ent, with photon indices of 3.47 + 0.07 and 2.80 + 0.04 respec- 
tively. We then investigated these differences through more de- 
tailed spectral modelling. From Table 1 one can see that the sim- 
ple powerlaw model gives a rather poor fit for observation 295. 
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Spectral PL-model fitting for GRS 1915+105 



1.2-12 keV ASM light-curve 




Fig. 3. Basic source characteristics : model photon index versus 
mean flux per observation. The observations that are discussed 
in the text are highlighted in red. 
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Fig. 4. ASM and SPI light-curves for observations 295 and 423. 
For both observations the source exhibits low flux and very little 
variabiHty in the ASM band. In the 20-50A:ey SPI band we 
observe similar stability on the science-window timescale with 
intermediate and high average flux respectively. 



More precisely, we noticed that above 100 keV all spectral points 
are located above the fitted powerlaw. Pure thermal comptoniza- 
tion models like COMPTT (where the comptonized spectrum is 
completely determined by the plasma temperature and its optical 
depth, Titarchuk 1994) can thus be ruled out (x^/v - 72/26, see 
table 2). Assuming that the low-energy part is nonetheless pro- 
duced through thermal comptonization, one needs to add a fur- 
ther spectral component. We thus chose to add a powerlaw as a 
phenomenological description of the observed high-energy tail. 
Given the error amplitude above 100 keV, we arbitrarily fixed the 
photon index to 2.0. The resulting fit (kTe - 16.3^q g keV and t - 
0.5^ is in very good agreement with our data (x^/v — 17/25) 
and the F-TEST indicates a probability of =s 10"^ that this im- 
provement has been a chance event. As a last step we applied 
the Poutanen & Svensson (1996) COMPPS model which de- 
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Fig. 5. ASM and SPI Ught-curves for observations 246 and 368. 
During the latter a correlation (on the alks science-window 
timescale) between the 1.2-12 and 20-50A;ey bands is no- 
ticeable, hinting that both bands are tracing the evolution of the 
comptonized emission. 



scribes comptonization from a hybrid thermal/non-thermal elec- 
tron plasma. This iterative scattering method assumes a pow- 
erlaw distribution (oc -y'^') up to a Lorentz factor y,„,„ below 
which the plasma thermalizes to a Maxwellian distribution. We 
fixed a spherical geometry and assumed the reflection compo- 
nent to be insignificant, enabling a straightforward comparison 
with COMPTT+PL. Note that COMPPS poses more constraints 
than COMPTT+PL, as it requires both components (thermal and 
non-thermal) to be linked, whereas the former does not. Equally 
good fitting results (x^/v - 17/25) for both models indicate that 
the observed 20-500fcey emission is most likely to originate 
from thermal and non-thermal comptonization processes. 

For observation 423, the basic powerlaw fit is clearly unac- 
ceptable (x^/v = 83/27) due to the marked curvature around 
50 keV, leaving clear evidence for thermal processes. Replacing 
the powerlaw by a pure thermal component improved the fit and 
gave;t'^/v = 42/26. However, this interpretation is not able to 
account for the observed emission above 200 keV, again suggest- 
ing the presence of an additional component. We thus keep the 
hybrid comptonization models (x^/v - 37/25 for COMPTT+PL 
andx^/y = 37/24 for COMPPS) as our preferred description for 
observation 423. 

All fitting results are summarized in Table 2 and discussed 
in section 4. 



3.2. Observations 246 and 368 

Both observations have already been intensively studied. 
Observation 246 is discussed by Rodriguez et al. (2008b) who 
conducted detailed spectro-temporal analysis by means of multi- 
wavelength coverage. They found that the source showed peri- 
odic X-ray cycles, identified as alternations of v and p classes 
(see Belloni et al. 2000 for definitions). INTEGRAL observa- 
tion 368 was again part of a large multi-wavelength campaign, 
this time involving the Suzaku satellite (Ueda et al. 2006). These 
authors identified the observed X-ray variability pattern to be the 
signature of a high flux transition from a soft class x to class 
0. We find the same variability pattern in X-rays as in the 20 - 
50 keV SPI band (Figure 5), indicating that both bands are prob- 
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ably sampling the temporal behavior of the same component. 
For observation 246 there seems to be a similar correspondence, 
although less clear due to lower variability amplitude in the SPI 
band. Both spectra are fitted with thermal + non-thermal comp- 
tonization models (COMPTT+PL and COMPPS) which provide 
the best agreement to the data (see Table 2). Compared to the 
joint JEMX/ISGRI spectrum from observation 246 (Rodriguez 
et al. 2008b), our SPI data are found to have almost the same 
spectrum than the one produced by these authors from only 
stable low/hard X-ray intervals (Interval I in Rodriguez et al. 
2008b), whereas the very short X-ray spikes originating from 
the disk (Intervals II and IV) are washed out. This shows that the 
SPI spectra of the one day averaged comptonized emission are 
not influenced by any peculiar outburst-behavior on very short 
time-scales. Likewise, a comparison with simultaneous Suzaku 
observations leads to similar conclusions for observation 368. 



3.3. Cross-comparison 

When comparing the results of these four high-lighted observa- 
tions, we notice that the X-ray variability and the spectral be- 
havior at higher energies (> 20 ^eV) seem to be completely un- 
correlated. Except for a slight difference in flux, the SPI spectra 
from observations 368 and 295 (illustrated in Figure 6) show 
the same characteristics, whereas the temporal behavior in X- 
rays is seen to be dramatically different (compare Figure 4 left 
and Figure 5 right). A similar high-energy correspondence can 
be observed between observations 246 and 423 (Figure 7), al- 
though again no similarities are found in X-rays. We thus see 
that for similar X-ray variability, the 20 - 500 A:e V spectra can be 
different and reciprocally that for similar ( ~ one day averaged) 
high-energy spectral behavior, the source may exhibit very un- 
like \.2-\2keV activity. 

3.4. Composite spectra 
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Fig. 6. COMPTT+PL best fits for observations 295 (red stars) 
and 368 (blue diamonds). Both observations have very similar 
spectral shapes. 



In order to further investigate the high-energy tail of GRS 
1915+105, we decided to group data from observations with 
similar spectra. Improved statistics should allow us to put bet- 
ter constraints on the parameters of the fitted models. Hence, 
we isolated two opposite groups in the spectral index versus 
flux diagram (Figure 3) for which we generated composite spec- 
tra. The first group is characterized by a rather low 20-50A:ey 
flux (^ 200 mCrab) and a very soft spectral shape (F ^ 3.45); 
hereafter we will call it the soft sample. The second group on 
the other hand has hard colors (F « 2.90) and high average 
flux (» 330 mCrab) in the 20-50A:ey band; it will accordingly 
be called the hard sample. We note that the previously pre- 
sented observations 295/368 and 246/423 are representative of 
each group, respectively. As some observations had intermediate 
characteristics, we excluded them from the regrouped spectra in 
order to avoid mixing different high-energy patterns. As a result, 
our two samples are likely to describe the boundary comptoniza- 
tion states between which the source seems to be continuously 
switching. 

We fitted both composite spectra shown in Figure 8 with sev- 
eral models and summarize the results in Table 2. Each time the 
thermal + non-thermal comptonization models gave the best fit 
to the data and demonstrated the need for an additional high- 
energy component (see F-TEST values in Table 2). The compos- 
ite spectral analysis thus confirms the previously outlined inter- 
pretation, which will be further discussed in the next section. 



4. Discussion 

We interpret the 20-500 /teV emission of GRS 1915+105 as a 
combination of thermal and non-thermal comptonization. Given 
the similarities of the observed temporal properties, the ASM 
band is likely to sample the same emission component as the 
20-50A:ey SPI band. This is in agreement with the results of 
Done, Wardziriski & Gierliriski (2004) who show that except for 
ultrashort disc-dominated X-ray spikes the accretion disc has no 
significant effect above 3 keV (see also Rodriguez et al. 2008b). 
The spectra obtained during observations 295 and 368 show that 
for different X-ray classes (i.e. different temporal behavior of 
the accretion flow) the high energy spectra can be very simi- 
lar. Conversely, observations 295 and 423 show that within sta- 
ble emission episodes, there can be significant differences in the 
spectral behavior of the source. Assuming that the 3-50keV 
emission originates from the comptonizing corona (see Malzac 
2007 or Done et al. 2007 for details), this clearly indicates that 
there is no correlation between the temporal variability and the 
macroscopic properties of the comptonizing flow. As they give 
better statistics above 100 ^eV, we will use our composite spec- 
tra to further discuss both comptonization components and put 
our findings in a more general context. 
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Table 1. Log of the SPI observations and basic source characteristics. Temporal behavior is parameterized in terms of average flux and fractional 
rms whereas the photon index T andx^/v are the parameter and quality of a simple powerlaw model fit. 
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Fig. 7. COMPTT+PL best fits for observations 246 (red stars) 
and 423 (blue diamonds). 



4.1. The soft sample 

From spectral modelling we find the intersection point of the two 
components to be around 60 keV. At higher energies, the non- 
thermal processes dominate. Both the thermal and non-thermal 
components have approximately the same luminosity above 
2QkeV. If these two components originate from the same pop- 
ulation (COMPPS results), the electrons are found to thermal- 
ize around Lorentz factor y,„,„ w 1 .30. The mildly hot Compton 
cloud (kT,, = 14.72 + 0.15 keV) is found to be marginally opti- 
cally thick (t - 2.7^y 2j) which is the expected configuration for 
black-hole binaries in the very high state (Done & Kubota 2006). 
We suggest that our soft composite spectrum gives a canonical 
high-energy representation of GRS 1915+105 in low coronal lu- 
minosity states. 

4.2. The hard sample 

Concerning the hard sample, the luminosity of the non-thermal 
component is found to be roughly the same (compatible within 
the uncertainties) as for the previously discussed soft sam- 
ple. However, the intersection point of the two components is 
now around llO^eV (which corresponds to a Lorentz factor 
•y„„„ ~ 1 .39), showing that the main difference lies in the prop- 
erties of the comptonizing thermal electrons of the corona. The 
plasma is found to be either hotter for similar optical depths or 
optically thicker for similar electron temperatures (or a mixture 
of both), thus enhancing the higher observed 20-50A:ey flux. 
Even though this situation cannot be completely resolved due to 
the observational kT -t degeneracy, spectral fits with COMPPS 
indicate that it is likely that there has been a significant in- 
crease in opacity, whereas electron temperature remains around 
15 - IQkeV. In any case, this does not affect the estimation of 
the total 20-500^ey luminosity issued from thermal Compton 
scattering, which is found to be enhanced by a factor of 2 + 0.04 
in comparison with the soft sample. We interpret our spectrum 
as a typical representation for the high-energy emission of GRS 
1915+105 in high coronal luminosity states. 
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Table 2. Spectral fitting results for the high-lighted observations. For each spectrum, the first row corresponds to a simple thermal comptonization 
model (COMPTT), the second row contains the COMPTT + PL parameters and the last gives the variables of the COMPPS model. The seed 
photon temperature is fixed to l.OfceV and 1.5keV for the harder and softer spectra respectively. COMPTT errors are calculated on kT^ while 
fixing T to its best fit value. For the COMPPS model, the geometry parameter is fixed to zero (spherical) and neither reflection nor ionization is 
taken into account. The error bars are simultaneously evaluated on kT^ and r while we fixed 7^i„ and to their best fit values. Fcomp,, denotes the 
integrated 20-500 keV energy flux from the thermal component and is given in units of xlO^'^ erg / cm^ / s. Kpi is the flux normalization at 100 feV 



of the non-thermal power law and is given in units of xlO ^photonscm ^keV 



5. Summary and Conclusions 

We have conducted detailed high-energy spectral analysis of the 
microquasar GRS 1915-H105 using all available SPI data from 
July 2004 to May 2006. We presented here our observational 
results on the source's high-energy output. We can summarize 
our findings as follows: 

We found that the ~ one day averaged 20-500 A:ey spectral 
emission is always between two boundary states, hard and soft, 
which we illustrated through spectral modelling. We confirm that 
in the INTEGRAL-SPI data we observe no high-energy cutoff 
for GRS 1915-1-105 (Fuchs et al. 2003, Rodriguez et al. 2008b). 
We suggest that the high-energy cutoff from thermal comp- 
tonization is drowned out by an additional non-thermal compo- 
nent. We found the non-thermal component to be statistically 
required in both composite samples. The spectral differences we 
observed in hard X-rays (20 - 50 keV) are most likely to be cou- 
pled to the evolution of the thermal electron plasma. The bolo- 
metric luminosity calculated from the thermal component varies 
by a factor of 2 between soft and hard samples. In contrast, the 
obtained fits indicate that the non-thermal component is fairly 
stable. This implies that both components are not necessarily 
linked, i.e. they could originate from dissociated electron pop- 
ulations. Hence as suggested by Rodriguez et al. (2008b), the 
non-thermal component might be from emission from the jet. 
Given the length of the high-energy observations (SPI » 3 days 
or OSSE 15 days), it is difficult to investigate the connection 
between the various X-ray classes and the high-energy spec- 
tra. We pointed out that there is no direct correlation between 



the observed X-ray variability patterns and our 20 - 500 feeV SPI 
spectra. This shows that the macroscopic properties of the comp- 
tonizing thermal electrons evolve independently from the tem- 
poral behavior of the source, i.e. independently from the fluctu- 
ations of the accretion flow. These aspects will be adressed in a 
subsequent paper using multiwavelength coverage. In summary, 
our results give a recent high-energy picture of GRS I9I5-HI05 
and once more underline the complex nature of the accretion 
processes operating around this archetypal microquasar. 
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Fig. 8. COMPTT+PL best fits for the composite spectra, plotted 
in red stars and blue diamonds for the soft and hard samples 
respectively. These spectra illustrate the boundary characteristics 
of the one day averaged high-energy emission of GRS 1 9 1 5 + 1 05 . 
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